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Abstract —A newly proposed chemical-reaction-inspired meta- 
heurisic, Chemical Reaction Optimization (CRO), has been 
applied to many optimization problems in both discrete and 
continnous domains. To alleviate the effort in tuning parameters, 
this paper reduces the number of optimization parameters in 
canonical CRO and develops an adaptive scheme to evolve them. 
Our proposed Adaptive CRO (ACRO) adapts better to different 
optimization problems. We perform simulations with ACRO on a 
widely-used benchmark of continuous problems. The simulation 
results show that ACRO has superior performance over canonical 
CRO. 

Index Terms —Chemical Reaction Optimization, continuous 
optimization, adaptive scheme, metaheuristic, evolutionary algo¬ 
rithm. 


1. Introduction 

Optimization techniques are frequently used in science 
and engineering research and development. For example, the 
problem of computing the contact force in mechanics can 
be regarded as a quadratic programming problem Q. In 
economics, firms are usually expected to maximize their profit, 
while customers aim to maximize their utilities. So the asset 
price can be modelled by optimization techniques i2. We 
also encounter numerous optimization problems in our daily 
lives, such as finding the shortest route to a destination while 
avoiding traffic congestion El, and assigning various tasks to 
different time slots to maximize the amount of work done 
before the deadline 0. 

Chemical Reaction Optimization (CRO), proposed in 0, 
is a chemical-reaction-inspired general-propose metaheuris¬ 
tic. CRO is variable population-based and it was initially 
designed to solve combinatorial problems. Then Lam et al. 
proposed a Real-Coded CRO (RCCRO) 0, which aims to 
solve continuous optimization problems. In this paper, with 
abuse of notation, we use CRO to refer to both CRO in 0 
and RCCRO in 0. In a chemical reaction, the reactants tend 
to release the excessive energy to the environment to approach 
relatively lower energy states. CRO utilizes this tendency 
and incorporates the chemical reaction ideas to construct an 
optimization algorithm. 

CRO has been applied to solve many optimization prob¬ 
lems. Lam and Li adopted CRO to solve Quadratic Assign¬ 
ment Problem (QAP), Resource-Constraint Project Scheduling 
Problem, and Channel Assignment Problem in 0. The first 
two are classical combinatorial problems while the last one 
is a practical NP-hard combinatorial problem according to 
0. They are used as benchmark problems for many other 
algorithms 0||9l and CRO outperforms many of such algo¬ 
rithms 0. Xu et al. ifTOll applied CRO to a task scheduling 


problem in grid computing, which is a multi-objective NP- 
hard optimization problem. Lam et al. addressed the cognitive 
radio spectrum allocation problem with CRO in ifTTIl . Yu et 
al. proposed a sensor deployment problem for air pollution 
detection and solved it with CRO m. All these problems 
mentioned are NP-hard combinatorial optimization problems 
and CRO was shown to be effective and efficient in solving 
these problems. Lam et al. also studied the convergence 
characteristics of CRO in ini. 

Lam et al. then proposed RCCRO to tackle continuous 
problems 0. They made three major modifications on the 
structures of molecule and algorithm to make CRO applicable 
to continuous problems. A classic benchmark set with 23 
continuous functions were adopted for testing and the simula¬ 
tion results indicated that RCCRO outperformed most of the 
other Evolutionary Algorithms (EAs) compared. Lam et al. 
proposed a population transition problem for peer-to-peer live 
streaming and applied CRO to solve the problem lfl4ll . Yu et al. 
proposed a CRO-based Artificial Neural Network (CROANN) 
algorithm to train neural networks (NN) for classification iflSll . 
CROANN also outperformed most previously proposed EA- 
based NN training algorithms. Yu et al. also studied different 
perturbation functions for RCCRO in ifThll . They compared 
the Gaussian distribution function used in 0 with three other 
distribution functions for the 23 benchmark functions and 
provided guidelines for further RCCRO research. 

According to our experience, proper parameter setting is key 
for CRO to solve optimization problems efficiently. However, 
we usually tune the parameters according to our experiences 
and by trial and error. Moreover, there are no theoretical results 
on parameter setting for EAs in general. To avoid the time- 
consuming parameter tuning process, we may pursue adaptive 
approaches. We developed a deterministic parameter control 
scheme for CRO in 0, but this scheme is still lacking as it 
only modifies the StepSize parameter of CRO deterministically 
and there is no feedback to inform the modification. On 
the contrary, adaptive parameter schemes have been proposed 
for some other EAs and they allow one to minimize the 
effort on parameter tuning and focus on the formulation of 
the optimization problem and the design of the algorithm 
operators. 

This paper is focused on developing an adaptive scheme for 
CRO. This new version of CRO with the adaptive scheme is 
called Adaptive CRO (ACRO). ACRO reduces the number of 
parameters defined in canonical CRO from eight to three, and 
makes them adaptive to feedback information in the course of 
searching. To evaluate its optimization performance, we use 
ACRO to solve a set of benchmark functions mainly selected 
from the latest IEEE CEC Real-Parameter Single Objective 
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Optimization Competition problem set ini. 

The rest of this paper is organized as follows. In Section 
im we briefly introduce the framework of CRO. We present 
ACRO with parameter adaptation in Section |III1 and Section 
HVl describes the benchmark problems used to examine the 
performance of ACRO and the simulation configuration. The 
comparisons with the canonical CRO are presented in Section 
[Vl We conclude the paper in Section |Vl] 

II. Chemical Reaction Optimization 

In this section, we introduce the general framework of CRO. 
We consider a number of molecules in a container, with an 
attached central energy buffer. 

A. Molecule 

A molecule is the basic operating agent of CRO. CRO 
controls and manipulates a set of molecules systematically to 
explore the solution space. Each molecule is characterized by 
attributes such as the molecular structure, potential energy, and 
kinetic energy. We relate chemical reactions to optimization 
by assigning mathematical meanings to these attributes. In 
a typical implementation of CRO, we include the following 
attributes: 

1) Molecular Structure: We denote the molecular structure 
by uj, which represents the solution carried by the molecule. 
The detailed structure of uj is problem-dependent, which 
means that the configuration of uj is usually determined by the 
encoding scheme of the problem. For different optimization 
problems, uj can be a vector, a matrix, a graph, or a string 

M- 

2) Potential Energy: We use to denote the potential 
energy (PE) of the molecule with molecular structure uj. In 
CRO, represents the objective function value of the 
solution carried by the molecule, i.e., f{uj). So our goal of 
optimization is to reduce as much as possible. Since 
“Every reacting system seeks to achieve a minimum of free en¬ 
ergy” mi, the chemical reaction process drives the molecules 
toward the lowest energy state on the potential energy surface 
(PES). 

3) Kinetic Energy: We use KE^^ to denote the kinetic 
energy (KE) of the molecule with molecular structure uj. In 
CRO, KE^ is regarded as the tolerance of the molecule to 
accept a worse uj' with higher PE^i than the existing PE^. 
The main purpose of introducing KE to CRO is to avoid the 
algorithm from getting stuck in local optima. This will be 
elaborated in Section Hl-CI 

B. Energy Conservation Law 

Energy transformation between different forms and conser¬ 
vation of energy are features which distinguish CRO from 
other metaheuristics 0. The law of conservation of energy 
states that although energy is allowed to transform between 
types, the total energy in an isolated system (i.e. the container 
in CRO) shall remain constant. We use E,otai and to 

denote the total energy in the container before and after an 
elementary reaction, respectively. Similarly, we denote the en¬ 
ergy in the central energy buffer before and after an elementary 


reaction by Etuffer and respectively. In an elementary 

reaction, consider that the original molecules cui, a; 2 , 
are transformed into new molecules wj, The total 

energy before the change is given by 

E,otal =Ehuffer + PE^j^ + KE^^^ + PE^^ -f KEi^^ -f ... 

+ PE^^ +KE^^, 

and it will remain constant, i.e., 

E'total =E'buffer + + PE^'^ -f EE^' -I-... 

+ PEoj'^ -f KEi^i^ (2) 

—Efotal • 

The transformation among Ehuger, PE, and KE for CRO is 
stated in 0 and interested readers can refer to m for more 
details. 

C. Elementary Reactions 

There are four types of elementary reactions defined in 
CRO, including on-wall ineffective collision, decomposition, 
inter-molecular ineffective collision, and synthesis. We classify 
the elementary reactions by the number of molecules involved 
(see Table IJl. For uni-molecular reactions, only one molecule 
is involved in the reaction as the reactant. For inter-molecular 
collisions, two molecules act as the reactants of the reaction. 
An elementary reaction may or may not change the total 
number of molecules in the system. Constant population size 
collisions do not change the number of molecules in the 
container, while variable population size collisions increase 
or decrease the number of molecules by one. A uni-molecular 
reaction is triggered when a molecule collides with a wall of 
the container, while an inter-molecular reaction corresponds 
to two molecules colliding with each other. Moreover, only 
successful reactions (subject to energy constraints) will modify 
the molecular structures of the involved molecules. CRO mod¬ 
ifies the molecular structures of the molecules to explore the 
solution space with these elementary reactions. The detailed 
information about these elementary reactions can be found in 
0 and d). 

III. ACRO 

In this section we propose our new CRO with parameter 
adaptation. 

Canonical CRO has eight parameters, namely, initial popula¬ 
tion size (iniPopSize), initial molecular kinetic energy (iniKE), 
initial central energy buffer (iniBuffer), molecular collision 
rate (CollRate), energy loss rate (LossRate), decomposition 
occurrence threshold (DecThres), synthesis occurrence thresh¬ 
old (SynThres), and perturbation step size (StepSize). These 
parameters cooperate to control the performance of the algo¬ 
rithm. 

In our proposed new CRO algorithm, we modify six 
parameters of CRO, including iniKE, iniBuffer, LossRate, 
DecThres, SynThres, and StepSize. We categorize them into 
three classes: energy-related, reaction-related, and real-coded- 
related as shown in Table HIl We will discuss our modifications 
in the following subsections. In general, we will simplify the 




TABLE I 

Classification of Elementary Reactions 



Uni-Molecular Collision 

Inter-Molecular Collision 

Constant Population Size Collision 

On-wall Ineffective Collision 

Inter-Molecular Ineffective Collision 

Variable Population Size Collision 

Decomposition 

Synthesis 


TABLE II 

CRO Parameter Categories 


Category 

Parameter 


iniKE 

Energy Related 

iniBuffer 

LossRate 

Reaction Related 

DecThres 

SynThres 

Real-Code Related 

StepSize 


design of iniKE, iniBujfer, and LossRate. Then we will use a 
new parameter to replace both DecThres and SynThres. The 
new parameter ChangeRate, along with a newly proposed el¬ 
ementary reaction selection scheme, will make the occurrence 
of elementary reactions adaptive. Finally we will employ a 
modified “1/5 success rule” to make StepSize adaptive. The 
remaining parameters {iniPopSize and CollRate) retain the 
original dehnitions in CRO presented in Q. 

A. Energy-Related Parameters 

Parameters iniKE, iniBujfer, and LossRate cooperate with 
each other to control the convergence of CRO lO. iniKE and 
iniBuffer determine the total energy in the system before the 
algorithm starts, while LossRate controls the rate of trans¬ 
forming KE to the central energy buffer in on-wall ineffective 
collisions. 

In CRO, KE controls the tolerance of the molecule to accept 
a worse structure a. So a large KE may potentially lead 
to a situation in which the molecule keeps accepting worse 
solution structures, resulting in bad optimization performance. 
Meanwhile, small KE may cause the molecule to get stuck 
in local optima since the molecule does not have enough 
tolerance to jump out of local optima. As stated in the 
energy of the central buffer and PE in molecules can be 
transformed into KE, and this process is partially controlled by 
LossRate in on-wall ineffective collision. So iniKE, iniBuffer, 
and LossRate also impact the performance, and proper values 
of these parameters are very important in controlling the 
convergence of the algorithm. 

1) iniKE: This parameter is used to control the initial KE of 
molecules. When the algorithm starts and the initial population 
is generated, all the molecules in the population are assigned 
with KE according to iniKE. The selection of iniKE shall be 
made based on the characteristics of the objective function 
value. 

In ACRO, we first generate a population of molecules 
popinit = -^miniPopSize} by choosing feasible solu¬ 

tions in the solution space randomly. Their corresponding PEs 
are determined. We label the molecular structures of the two 
molecules with the largest and the smallest PE by uji and ujg in 


popi„i,, respectively. Then we define the initial kinetic energy 
of all the molecules as 

iniKE = (PEcj, — PE^i^ ) x iniPopSize (3) 

and this initial KE rule is applied whenever a new molecule 
is generated in subsequent decompositions. 

2) iniBuffer: The central energy buffer Ebuffer accumulates 
energy in on-wall ineffective collisions, and energy is with¬ 
drawn in decomposition HD. As Ebuffer is initiated according to 
iniBuffer, this value is very important in determining whether 
a decomposition is successful, especially in the early stage of 
the search. 

Decomposition is often manipulated with operators per¬ 
forming vigorous modihcations to the molecular structures, 
assisting molecules to jump out of local optima. But a de¬ 
composition that occurs too early may potentially hamper 
the convergence of the algorithm, especially in those cases 
in which the decomposition modifies the molecules that are 
still evolving, i.e. not stuck in a local optimum. This hinders 
ACRO from performing an effective optimization search as 
the molecule will switch to another region of PES instead of 
performing a local search. In other words, this may lead to 
missing a possible global optimum in a particular region. As 
stated previously, the success rate of decomposition is partially 
controlled by Ebuffer- We can reduce the initial energy in the 
central buffer to suppress the decomposition at the early stage 
since a small central buffer, initialized with a small iniBuffer, 
renders it less likely for decomposition to occur ||6l. When 
the molecules have collided and explored the solution space 
for some time, Ebuffer is then gradually increased due to on- 
wall ineffective collisions, and the molecules are more likely 
to get stuck in the local optima. Since Ebuffer is getting larger, 
more decomposition will be successful and the molecules can 
explore other parts of the solution space without getting stuck. 

Therefore, according to the previous analysis, we assign 
zero to iniBuffer to suppress decomposition in the early stage 
of searching in ACRO. The idea of setting the initial buffer 
to zero was also supported numerically by the simulations in 
EOl and i). 














3) LossRate: Parameter LossRate is used in the energy 
transformation process of the on-wall ineffective collision. 
When a molecule is involved in a successful on-wall inef¬ 
fective collision, a portion of the energy that the molecule 
originally holds will be transferred to Ebuffer- 

In a successful on-wall ineffective collision, the excessive 
energy E^xcess is divided into two parts. One part becomes the 
KE of the new molecule according to 

KE^/ = Egxcess ^ ^ ( 4 ) 

where w' is the new molecular structure, and ^ is a random 
number generated from the interval [LossRate, 1]. The remain¬ 
ing energy E^xcess — KE^,' is transferred to Ebuffer- Hence, Loss- 
Rate affects the KE a molecule possesses after a successful 
on-wall ineffective collision. 

In many previous works of CRO, LossRate is set to be 
a very small number (e.g., 0.1 in lb), 0.2 in ifT^ I. In our 
ACRO algorithm, the parameter LossRate is designed to be 
a molecular attribute, instead of the original global attribute. 
The value for LossRate of each molecule in the system is 
approximated by a modified folded normal distribution. This 
distribution is generated from a normal distribution with a 
mean value of 0 and standard deviation of 0.3. All negative 
values in this normal distribution is “folded” over by taking 
their absolute values, and all values larger than 1 are consid¬ 
ered as 1. By doing so, we attempt to make the energy loss 
rate a small enough number, but without losing the original 
searching characteristics provided by large LossRate. 

B. Reaction-Related Parameters 

In CRO, we use CollRate, DecThres, and SynThres to 
control the ratio among the four kinds of elementary reactions. 
We first apply CollRate to determine whether a uni-molecular 
collision or an inter-molecular collision would happen in the 
current iteration. Then one or two molecules are randomly 
chosen from the population as the input molecule(s) to an 
elementary reaction. For a uni-molecular collision, we com¬ 
pare the “inactive degree” 0 of the molecule with DecThres. 
If the inactive degree is larger than DecThres, decomposition 
will occur. For an inter-molecular collision, the KEs of the two 
input molecules are compared with SynThres. If both of the 
KEs are smaller than SynThres, synthesis will occur. However, 
this scheme may result in an imbalance between the number 
of decompositions and syntheses when inappropriate values 
are set to CollRate, DecThres, and SynThres. For example, if 
we have too many syntheses, the population size is likely to 
be reduced to one most of the time during simulation. Then it 
is highly likely for the molecules to get stuck in local optima 
and the population diversity is no longer maintained. On the 
other hand, if there are many decompositions, too many new 
molecules will be generated. The limited number of function 
evaluations will be exhausted without thoroughly searching the 
potential regions. Then, CRO may turn into a random search. 
For simplicity but without loss of generality, we assume that 
the population diversity can be maintained when the number of 
molecules is stable in the course of searching, with population 
size similar to the initial value, i.e., iniPopSize. 


In ACRO, we introduce a new parameter ChangeRate to 
replace the original parameters DecThres and SynThres to con¬ 
trol the frequency of decompositions and syntheses. Generally 
ChangeRate describes the probability of having decomposition 
and synthesis. For example, if we set ChangeRate = 0.01, 
then we will roughly have 1000 x 0.01 = 10 decompositions 
or syntheses among 1000 elementary reactions. However, 
ChangeRate alone is not enough to control the population 
size and maintain the population diversity. As in the previous 
example, it is possible that all 10 reactions are syntheses. Then 
if the iniPopSize is set to 11, the population size will finally 
be reduced to one. In such a case, the population diversity 
is still not maintained. In order to handle this problem, we 
shall increase the occurrence rate of decomposition when the 
population size is smaller than IniPopSize, and increase the 
rate of synthesis when the population size is large. To do this, 
we introduce two population feedbacksand/^^^ to control 
the occurrences of decomposition and synthesis, respectively. 

In ACRO, in each iteration, the algorithm will first decide 
whether it is a variable population-size collision or constant 
population-size collision using ChangeRate. Then the algo¬ 
rithm will decide whether a decomposition or synthesis shall 
happen. We define the population feedback term as 

curPopSize — iniPopSize 
iniPopSize 

Then the decomposition and synthesis feedbacks are defined 
as ^ 

fdec = 2 ^ 

and ^ 

fsyn = 1 -fdec = 3 ^ +fpop)^ (7) 

respectively, where curPopSize is the current population size. 

and determine the probability the current iteration is 
a decomposition and a synthesis, respectively. 

The population feedback term^^^ is positive when curPop¬ 
Size is larger than iniPopSize. A positive will encourage 
the occurrence of synthesis by increasing its probability and 
decreasing the probability of decomposition, and vice versa. 
If the current population size equals iniPopSize, = 0 
and = 0.5, which means that decomposition and 

synthesis are equally likely to occur. If the current population 
size is one and iniPopSize is not one, = —1 and the 
probability of having synthesis is zero. This matches the 
requirement of synthesis that at least two molecules in the 
population are needed. If iniPopSize and curPopSize are both 
one, despite= 0.5, we stipulate that only decomposition 
can occur. To conclude,cooperates with/^^^ to control the 
population size and keeps its fluctuations within a small range 
by controlling the relative occurrences of decomposition and 
synthesis. 

C. Real-Code-Related Parameter 

For continuous optimization problems, there are infinite fea¬ 
sible values for each variable. In CRO, we apply the zero-mean 
Gaussian distribution in the neighborhood searching operator 
for continuous optimization problems. Parameter StepSize is 



employed to control the variance of the Gaussian distribution, 
which determines “how far” the molecule can go in the on-wall 
and inter-molecular ineffective collisions E). This parameter 
is set to be a constant in CRO, and decreased at a predefined 
rate in the adaptive scheme Q. 

In CRO, the neighborhood search operator works as fol¬ 
lows. Given a molecular structure, we first randomly pick an 
element from it. Then we randomly generate a new number 
SA ^ -^(0, StepSize) according to the Gaussian distribution 
and add this number to the value of the element. Suppose the 
selected element is the i-th element in the solution, and its 
original value is Si. Then the new value s/ of this element 
after applying the neighborhood search operator is 

Si = Si -f SA- (8) 

In CRO, a proper value of StepSize is critical to the 
performance of the algorithm. While an overly large StepSize 
prevents the algorithm from converging, an unreasonably small 
StepSize will result in the algorithm getting stuck in local 
optima. The modification in StepSize is twofold; initial value 
assignment and its subsequent update. 

1) Initial Value of StepSize: In ACRO, we set an initial 
value to the StepSize and then update this parameter in the 
course of searching. For each element in a solution, there is 
an upper bound and lower bound to confine the feasible range. 
We set the initial StepSize for this element to be 

StepSize= {Boundupper,i - Boundiower,i)/‘2 (9) 

where StepSize ,■ is the initial StepSize for the i-th element 
in the solution, Boundupper,i and Boundiower are the upper and 
lower bounds for the i-th element. For simplicity, we omit i 
in the subscript of StepSize in the sequel. 

2) Evolution of StepSize: We adapt the “1/5 success rule” 
II 20 I to modify StepSize in the course of searching. This rule 
was originally stated as follows; 

After every n mutations, check how many successes have 
occurred over the preceding lOn mutations. If the number is 
less than 2n, multiply the step lengths by the factor 0.85; divide 
them by 0.85 if more than 2n successes occurred. 

The “1/5 success rule” has been applied to EP, ES and 
other EAs successfully ED and demonstrated outstanding 
performance. So we adopt this rule to manipulate the “step- 
size” parameter adaptively. However it cannot be employed 
directly in CRO since there is no “mutation” defined in CRO. 
In CRO, whenever an elementary reaction is successfully 
performed, new solution(s) are generated, and we have an 
“update” process at the end of each iteration 0 . In “update”, 
the algorithm compares the newly explored solutions with the 
best ever solution. If one of the new solutions is even better, the 
algorithm stores its PE (objective function value) and regards 
this solution as the best ever solution. In ACRO, we apply a 
modified “1/5 success rule” in “update”. 

The modification is as follows. The frequency at which 
StepSize is changed is kept the same, i.e., lOn updates, and the 
adaptation factor of StepSize is also kept constant at 0.85. n is 
set to be one hundredth of the maximum allowable function 
evaluation number (FELimit). The new update mechanism 
with “1/5 success rule” is presented in Algorithm [D 


Algorithm 1 ONUPDATENEWSOLUTION(w„eu;, W&est) 

1 : Compare the new solution uj„ew with the stored best ever 
solution LCbest- 
2 : if ujnew is better then 

3 ; ^best — ^new 

4 : This update is successful. 

5 : end if 

6 : for Every n updates do 

7 : if Previous lOn updates have more than 2n successful 

ones then 

8 : StepSize = StepSize / 0.81) 

9 : else 

10 : StepSize = StepSize x 0.85 

11 : end if 

12 : end for 


D. Summary 

In ACRO, we make StepSize adaptive. We remove Loss- 
Rate, DecThres, and SynThres and introduce a new parameter 
ChangeRate to control the occurrence of the elementary reac¬ 
tions. In addition, we formulate rules to assign initial values 
for iniKE and iniBuffer. So we reduce the total number of 
parameters from eight to three as indicated in Table [Till The 
decrease in the total number of parameters can reduce the 
effort to tune the parameter values. We will show that ACRO 
performs equally well or even better than the canonical CRO 
in the next section. 

IV. Experiment Setting 

In order to compare the performance improvement of our 
proposed ACRO algorithm over the canonical CRO in solv¬ 
ing global optimization problems, we conduct a series of 
simulations on a full set of different benchmark functions. 
As revealed in ll22l . many benchmark functions commonly 
employed for evaluating the optimization performance of 
evolutionary algorithms suffer from two major problems; the 
global optimal points are located at the center of the search 
space, and they are positioned along the coordinate axes, 
i.e., there is no coiTelation among different dimensions. We 
can shift and rotate the conventional benchmark functions 
to resolve these problems. We make a comprehensive test 
suite by employing 26 benchmark functions mainly selected 
from the latest CEC 2013 Real-Parameter Single Objective 
Optimization Competition El as well as the benchmark 
functions employed to test the performance of canonical CRO 
in 0 . All of these benchmark functions except /13 and /14 are 
shifted and a part of them are further rotated. The benchmark 
functions are listed in Table |IV] All benchmark functions are 
simulated on 30 dimensions. The search space and population 
initialization range are both defined to be [—100,100]. The 
global optimum value of these benchmark functions are all 0 , 
and all simulation results smaller than 10 “® are considered as 

0 El- 

Both ACRO and canonical CRO are implemented in Python 
2.7 on Microsoft Windows 7. All simulations are performed 
on a computer with an Intel Core 17-3770 @ 3.4GHz CPU. 







TABLE III 

Parameter Modifications in ACRO 


Parameter 

Initial Design in CRO 

Modification in ACRO 

MKE 


The initial KE is assigned according to the PE of the initial population. 
The formula is iniKE = (PEoji — PE^j^ ) x iniPopSize. 

iniBujfer 


The initial central energy buffer is set to zero. 

LossRate 

User assigns a proper initial value. 

This parameter is unique for each molecule, and the value is generated 
from a folded normal distribution. 

DecThres 

SynThres 


These t\vo parameters are replaced with a new parameter ChangeRate, 
to control the occurrences of decomposition and synthesis. 

StepSize 


I’he initial StepSize is set according to the feasible solution bounds. It 
evolves according to a modified “1/5 success rule”. 

iniPopSize 

CoURate 

No modification 


TABLE IV 

Benchmark Functions 


Solution Transformation 


Name 


E d 2 

i = l ^j = l 

max{\zi\, l<i<D} 
max{\zi\, 1 <i< D} 

+ - 1 )^) 

Y:,=\^(^00(z^ - ., + 1)2 + (., - 1)2) 


fl 

f2 

/3 

/4 

/5 

/6 

fl 

fs 

f9 

flO 

fll 

fl2 

fl3 

fl4 

/l5 

fl6 

fll 

fl8 

fl9 

f2Q 

f21^ 

f22^ 

/23^ 

/24^ 


-20 exp(-0.2^ — ^ “ exp(— 

+ 20 + e 


-20exp(-0.2i/ — 1 ) “ exp(— cos27r2,) 

+ 20 + e 

418.9829D - sin y^l z, I 

418.9829D - Yl^—i z, sin Y^l z, I 
- 10cos27rz, + 10) 

~ 10cos27rz, + 10) 

^^ = 1 4000 /Ji i 

n z'^ D 

^^ = 1 4000 ,i\ i 

sin2(7ri/i) + ” 1)^(1 + lOCsin^ y^_^^))] + 

{yn - 1)^(1 +sin2(27ryn)), Vi = 1 + -(^i + 1) 
sin2(7ri/i) + X^^_i [(Vi - 1)^(1 + 10(sin2 t/, + i))] + 

(l/n - 1)^(1 +sin2(27rT/n)), Wi = 1 + + 1) 

1 4 

- [sin2(37r2i) + “ 1)^(1 + sin^ (Sttz^+i )) + 

10 ^z-1 ^ 

(zn - 1)^(1 + ain2(27r2n,))] + H^^i - 5, 100, 4) 

— [sin2(37rzi) + “ 1)^ C1 + sin^ (37rz,+1) ) + 

10 ^i-l ^ 

(zn - 1)^(1 + sin2(27rzn))] +13^=1 100, 4) 

— [10sin2{7ri^l) + (vi - 1)^(1 + 10 sin^ {7ry,+1)) + 

(Vn - if] + Z;^=i ^(^i. 10, 100, 4), y, = 1 + ^(z, + 1) 

— [10sin2(7ryi) + “ 1)^(1 + 10 sin^ (7ry,+ i)) + 

(Vn - if] + ^(^i. 10, 100, 4), y, = 1 + -(z, + 1) 


5 = lVI(x - o) 

= = IVI(x - o) 

! = (x - o) X 0.1 
5 = M(x - o) X 0.1 
5 = (x - o) X 0.3 
5 = M(x - o) X 0.3 


5 = (x - o) X 0.32 

5 = M(x - o) X 0.32 

5 = X X 5 

5 = Mx X 5 

5 = (x - o) X 0.0512 


Shifted Sphere Function 
Shifted Schwefel’s Problem 1.2 
Shifted Rotated SchwefeTs Problem 1.2 
Shifted Schwefel’s Problem 2.21 
Shifted Rotated SchwefeTs Problem 2.21 
Shifted SchwefeTs Problem 2.22 
Shifted Rotated SchwefeTs Problem 2.22 
Shifted Rosenbrock's Function 
Shifted Rotated Rosenbrock’s Function 
Shifted Discus Function 

Shifted Ackley's Function 

Shifted Rotated Ackley's Function 

SchwefeTs Problem 2.26 
Rotated SchwefeTs Problem 2.26 
Shifted Ra.strigin's Function 


5 = ]VI(x — o) X 0.0512 Shifted Rotated Rastrigin’s Function 
5 = (x — o) X 6 Shifted Griewank's Function 


! = M(x - o) X 6 

5 = (x - o) X 0.1 


Shifted Rotated Griewank’s Function 


Shifted Levy’s Function 


: IVI(x — o) X 0.1 Shifted Rotated Levy's Function 


5 = (x — o) X 0.5 Shifted Penalized Function 1 


5 = M(x — o) X 0.5 Shifted Rotated Penalized Function 1 


: (x - o) X 0.5 


: M(x - o) X 0.5 


Shifted Penalized Function 2 


Shifted Rotated Penalized Function 2 


o is a shifting vector and M is a 

transformation matrix, o and AT can be obtained from 1171. 

1 

[ k(x — a)™ 

for X > a 

u{x^ a, /c, m) = < 

0 

for — a < X < a. 

1 


for X < —a 





















In order to reduce statistical errors and generate statistically 
significant results, each benchmark function is repeated for 51 
independent runs for each algorithm, which also satisfies the 
requirement of ini. In each run, we use the maximum number 
of function evaluations (maxFEs) as the termination criteria: 
100 000 maxFEs are employed for 10 dimensional tests and 
300 000 maxEEs for 30 dimensional tests ifTTl . 

To evaluate the performance improvement of our proposed 
ACRO algorithm over the canonical CRO, we compare the 
simulation results between these two algorithms. Lam et al. 
proposed four different CRO in 0, distinguished by dif¬ 
ferent constraint handling schemes, synthesis operators, and 
deterministic stepsize adaptation schemes. We use CRO/BP, 
CRO/HP, CRO/BB, and CRO/D to refer to RCCROl, RC- 
CR02, RCCR03, and RCCR04 in 0, respectively. Similarly, 
we consider the two different constraint handling schemes and 
two synthesis operators discussed in ||6|, and propose three 
variants of ACRO for performance comparison: 

1) ACRO/BP: This ACRO variant employs the boundary 
constraint handling rule (9) discussed in 0, and the 
“probablistic select” synthesis operator. 

2) ACRO/HP: This ACRO variant is based on ACRO/BP 
and the boundary constraint handling rule is the hybrid 
handling rule (10) in Q. 

3) ACRO/BB: This ACRO variant is based on ACRO/BP 
and the synthesis operator is replaced by BLX-0.5 in 
0 . 

The parameters for CRO are set according to the recommen¬ 
dation of 161 for solving multimodal problems, i.e. population 
size is 20, stepsize is 1, initial energy buffer is 10®, initial 
kinetic energy for molecules is 10^, molecular collision rate is 
0.2, kinetic energy loss rate is 0.1, decomposition threshold is 
1.5x10®, and synthesis threshold is 10. The stepsize adaptation 
interval for CRO/D is set to 100 and the change rate is 0.99. 
As we only have three parameters for ACRO, we set the 
population size and the molecular collision rate the same with 
CRO configurations, and ChangeRate as 10“^. 

V. Simulation Results 

In this section we present the simulation results of different 
variants of our proposed ACRO and canonical CRO described 
in Section|IVl We present the raw simulation result data as well 
as the statistical analysis based on the data. We also analyze 
the convergence of the different algorithms. 

The mean value obtained by all ACRO and canonical CRO 
algorithms on different benchmark functions are presented in 
Table [V] The mean values in bold font indicate superiority 
of the corresponding algorithm over the benchmark function. 
The simulation results indicate that ACRO generally gives 
superior performance over the canonical CRO algorithms. 
ACRO algorithms generate best results in 14 out of 24 30- 
D functions. 

From the simulation results in Table |V] the following key 
points can be observed: 

1) ACRO algorithms generally have better performance 
than CRO algorithms in all benchmark functions when 


comparing the mean simulation results, especially in uni- 
modal functions (fi-fio). Meanwhile, the superiority 
of ACRO in solving multi-modal functions is still very 
significant. While the deterministic stepSize adaptation 
scheme proposed in 0 has better overall performance 
than other CRO variants, ACRO algorithms outperform 
it by a large margin. 

2) While CRO has no significant preference on the bound¬ 
ary handling scheme selection, ACRO favors the bound¬ 
ary constraint handling rule (9) discussed in 0 more 
than the hybrid handling rule discussed in the same 
paper. This phenomenon can be observed in the bench¬ 
mark functions where ACRO/BP performs better than 
ACRO/HP, e.g. fr and /13. 

3) ACRO has no significant preference between the two 
different types of synthesis employed by ACRO/BB and 
ACRO/BP. These two algorithms perform similarly in 
all the benchmark functions tested. 

In terms of the computational time consumed by the pro¬ 
posed adaptation scheme, we observed that the additional time 
needed for ACRO algorithms is less than 5% of the computa¬ 
tion time of canonical CRO algorithms. As the computational 
overhead of the evolutionary algorithms on a given number 
of function evaluations with the same parameter setting is 
approximately constant when solving different optimization 
problems, we believe that the extra time used by the adaptation 
scheme can be considered negligible when solving problems 
with relatively long fitness evaluation time. 


VI. Conclusion 

In this paper, we proposed a new CRO-based metaheuristic 
called ACRO. To alleviate the parameter tuning process, we 
reduce the number of parameters defined in CRO. We also 
design an adaptive scheme for CRO. We manipulate six of 
the existing CRO parameters and successfully reduce the total 
number of parameters from eight to three. The advantages of 
ACRO over canonical CRO include: 

1) ACRO achieves superior optimization performance over 
canonical CRO within a similar amount of time. 

2) ACRO convergences faster than canonical CRO, and can 
effectively prevent pre-mature convergence. 

3) ACRO can significantly reduce the time needed to find 
a suitable combination of optimization parameters. 

These advantages can be observed in the simulation we con¬ 
ducted over a set of 24 benchmark functions. The simulation 
results demonstrate the superiority of ACRO. 

In the future, we will test the performance of ACRO on other 
practical problems. We will refine the settings for the proposed 
adaptive schemes. We can also include ACRO in an ensemble- 
based algorithm for real-parameter optimization Il23]lll24l . De¬ 
veloping ACRO-hybrid algorithms with other swarm intelli¬ 
gence algorithms, such as Particle Swarm Optimization 1251 
and Social Spider Algorithm l2^ . is also a potential research 
topic. 


TABLE V 

The Mean Simulation Results oe ACRO and Canonical CRO 


Problem 

ACRO/BP 

ACRO/HP 

ACRO/BB 

CRO/BP 

CRO/HP 

CRO/BB 

CRO/D 

/i 

O.OOOOe+00 

O.OOOOe+00 

O.OOOOe+00 

2.7374e-06 

2.6684e-06 

2.6009e4)6 

2.9454e-05 

/2 

O.OOOOe+00 

O.OOOOe+00 

O.OOOOe+00 

1.8280e-05 

1.9248e-05 

1.8690e-05 

1.4972e-04 

h 

1.9827e-06 

5.865 le-06 

1.4231e-06 

4.3676e-04 

5.3923e-04 

4.1359e-04 

1.4130e-01 

h 

6.0942e-01 

7.3579e-01 

I.0329e+00 

4.7625e-03 

4.7484e-03 

4.4301e-03 

1.2716e-02 

fd 

4.8004e+00 

7.1049e+00 

6.1221e+00 

5.8131e+01 

5.8909e+01 

5.3928e+01 

3.5403e+01 

h 

2.4722e-07 

2.1478e-07 

9.0862e-05 

5.2846e-04 

5.9767e-04 

6.2956e-04 

1.7779e-03 

fl 

8.4892e+01 

1.0700e+02 

8.7758e+01 

1.2944e+07 

6.8512e+06 

4.6210e+06 

2.2111e+01 

h 

2.4860e+01 

2.9761e+01 

2.4685e+01 

8.4063e+01 

8.4381e+01 

7.1735e+01 

6.8128e+01 

h 

2.2388e+01 

2.3080e+01 

2.2836e+01 

3.1328e+01 

2.9001e+01 

2.8922e+01 

3.0114e+01 

/lO 

1.7196e-07 

3.7136e-01 

4.2115e-06 

5.1422e-03 

5.6744e-03 

1.1732e-02 

1.9617e-01 

/ll 

8.4478e-02 

9.2534e-02 

7.8069e-02 

1.1503e+01 

1.1305e+01 

1.1177e+01 

2.1326e-04 

/l2 

2.3127e+00 

2.4233e+00 

2.3163e+00 

1.1388e+01 

1.1562e+01 

1.1231e+01 

3.0470e+00 

/l3 

2.9281e+02 

6.8552e+02 

2.3939e+02 

5.8231e+03 

6.0672e+03 

6.1206e+03 

O.OOOOe+00 

/l4 

4.1386e+03 

3.6792e+03 

3.8300e+03 

5.7527e+03 

5.5973e+03 

5.6308e+03 

8.6744e+02 

/l5 

6.9226e+00 

6.1607e+00 

4.8700e+00 

4.2171e+02 

4.3384e+02 

3.894 le+02 

1.4585e-04 

/l6 

1.6828e+02 

1.9685e+02 

1.6156e+02 

4.3566e+02 

4.5059e+02 

4.1632e+02 

2.4930e+02 

/l7 

3.0256e-03 

2.3141e-03 

1.5916e-03 

2.3064e+00 

2.1902e+00 

2.6168e+00 

1.0624e-02 

/l8 

5.8008e-04 

1.9328e-04 

2.9004e-04 

4.7648e-03 

1.9807e-03 

4.0542e-03 

9.1952e-03 

/l9 

O.OOOOe+00 

O.OOOOe+00 

O.OOOOe+00 

1.9428e+01 

1.2036e+01 

1.0251e+01 

1.8466e-08 

f20 

3.9330e+01 

4.308 le+01 

4.0080e+01 

3.7536e+01 

3.8191e+01 

6.9319e+01 

7.2789e+01 

/ 2 I 

4.0374e-08 

3.4969e-06 

9.4729e-10 

2.6110e-07 

2.3093e-07 

2.1121e-07 

1.5477e-06 

/22 

2.2444e+00 

2.2262e+00 

2.5946e+00 

2.9582e+01 

2.7701e+01 

2.8518e+01 

1.9143e+00 

/23 

1.8949e-03 

1.7269e-05 

1.2809e-03 

1.6633e+01 

1.5481e+01 

6.9264e+00 

1.5812e-07 

/24 

3.5586e+00 

3.5509e+00 

3.4288e+00 

2.1111e+00 

2.4530e+00 

2.5946e+00 

3.1687e+00 
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